Gravitational Shock Waves and Their Scattering in Brane-Induced Gravity 
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In this paper, following hep-th/0501028 we present a detailed derivation and discussion of the 
exact gravitational field solutions for a relativistic particle localized on a tensional brane in brane- 
induced gravity. Our derivation yields the metrics for both the normal branch and the self-inflating 
branch Dvali, Gabadadze and Porrati (DGP) braneworlds. They generalize the 4D gravitational 
shock waves in de Sitter space, and so we compare them to the corresponding 4D General Relativity 
solution and to the case when gravity resides only in the 5D bulk, and there are no brane-localized 
graviton kinetic terms. We write down the solutions in terms of two-variable hypergeometric func- 
tions and find that at short distances the shock wave profiles look exactly the same as in 4D 
Minkowski space, thus recovering the limit one expects if gravity is to be mediated by a metastable, 
but long-lived, bulk resonance. The corrections far from the source differ from the long distance 
corrections in 4D de Sitter space, coming in with odd powers of the distance. We discuss in detail 
the limiting case on the self-inflating branch when gravity is modified exactly at de Sitter radius, 
and energy can be lost into the bulk by resonance-like processes. Finally, we consider Planckian 
scattering on the brane, and find that for a sufficiently small impact parameter it is approximated 
very closely by the usual 4D description. 
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I. INTRODUCTION 



Among many recent braneworld proposals, a partic- 
ularly interesting idea called brane-induced gravity was 
put forth by Dvali, Gabadadze and Porrati (DGP) Q. In 
this approach the gravitational force between masses on 
the brane should arise from the exchange of an unstable 
bulk graviton resonance, which falls apart at very large 
scales, instead of a stable zero-mode graviton governed 
by the usual four-dimensional (4D) General Relativity 
(GR). Far from sources the law of gravity is modified, 
and changes from the 4D to a fully higher-dimensional 
one. This is similar to the approach of 2|, which is 
known to suffer from problems with ghosts Q (see how- 
ever for a possible cure of the ghost-induced mal- 
adies). The brane-induced gravity setup of DGP may 
be free of ghosts in perturbation theory [lj, and so it 
is an interesting arena to explore modifications of grav- 
ity at large distances. In naive perturbation theory, the 
modifications appear already at the level of linearized 
theory, where one finds that the tensor structure of the 
graviton-matter coupling changes from its form in 4D 
GR [l|. Indeed, around flat space, a small perturbation 



h^ u is sourced 



by - a-q^T \ , where a = 1/3 in- 



stead of 1/2, because of the extra scalar polarization of 
the massive graviton. Since this coupling does not de- 
pend on the graviton mass, the factor 1/3 would appear 
to persist even as m g — ► 0. This is an example of van 
Dam-Veltman-Zakharov (vDVZ) discontinuity 5], which 
could modify the predictions for the Solar system tests 
of general relativity. However, Vainshtein has argued 
that the naive perturbation theory for massive gravity 
breaks down in the limit m g — ► because it is a double 
expansion in both Gn and l/m g . Instead, the argument 
goes, one needs to expand around the curved background 



created by the source mass, where the curvature gener- 
ated by the source "screens" the scalar graviton, so that 
the modified perturbation theory smoothly reduces to 
conventional GR as m g — > 0. Basically, it is necessary 
to take into account the non-linear nature of gravity to 
properly restore gauge invariance which is required to 
ensure the decoupling of the massive modes as m g — > 0. 
A similar behavior has also been noticed in non-abelian 
gauge theories @- 

The issue of strong coupling in gravity with IR modifi- 
cations, and in particular its implications for DGP mod- 
els have been studied in pi la 111 liii liil \1M • The works 
El El 0] specifically address the longitudinal gravi- 
ton and its matter coupling, and suggest that brane ex- 
trinsic curvature may indeed play the role of a coupling 
controller. When a source mass is placed on the brane, 
it bends the brane. The brane extrinsic curvature which 
locally responds to the source may in turn help to tame 
the perturbation theory |12|. This interplay between 
the source-induced background curvature and the scalar 
graviton coupling may thus be a key ingredient of the 
"gravity filter" of It is therefore essential to explore 
in detail this interplay of the curvature and the coupling 
strength beyond perturbation theory in order to clarify 
the status of the effective 4D theory. However, the DGP 
equations are very difficult to solve exactly for compact 
sources localized on the brane 0, E| • 

In the previous paper E^] we have presented the first 
explicit example of an exact solution for a localized par- 
ticle in DGP models. We have obtained the gravita- 
tional field for a relativistic particle on a tensional brane 
in 5D, whic h g eneralizes 4D gravitational shock waves 
El El El El El El in de Sitter space El 0. 
In this paper we follow up with a detailed derivation of 
the shock wave metrics for both the normal branch and 
the self-inflating branch variants of DGP braneworlds 
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[25|. Focusing on the extreme relativistic limit, we find 
that the brane-localized terms are indeed a very efficient 
"gravity filter" . At short distances , th e wave profile be- 
haves exactly like in 4D GR 0, The deviations 
appear only far from the source, where they come as odd 
powers of the distance from the source in the units of the 
cosmological horizon length, as opposed to the 4D GR 
corrections in de Sitter space, which only come as even 
powers. We also look at gravitational scattering of rel- 
ativistic particles in DGP braneworlds, and find that to 
the leading order it behaves the same as in the 4D case. 

We note that our solutions do not provide a direct an- 
swer to the conundrum of strong couplings of the scalar 
graviton. However, we do find that the scalar graviton 
field is not turned on in these metrics. This would not 
come as a surprise in a weakly-coupled perturbation the- 
ory, where the scalar graviton will not be sourced by 
ultra-relativistic particles since to the leading order their 
source vanishes, T 7 ^ , = 0. In other words, in the rela- 
tivistic limit, one finds that the conformal symmetry is 
restored, and the solution is form-invariant under boosts. 
Hence the gravitational shock wave of 4D GR is also a 
solution in Brans-Dicke theory. In a more relevant case 
of very fast particles with non-zero rest mass, the cor- 
rections from the mass to the relativistic wave profile 
should be suppressed by the source mass-to-momentum 
ratio m/p <C 1. Our results show that decoupling of the 
scalar graviton persists in the relativistic limit in DGP, 
and so the dangerous strongly coupled mode identified 
in perturbation theory does not destroy the solutions at 
large momenta after all. This suggests that even for the 
particles with nonzero rest mass the gravitational field 
may be under control. Indeed, since local physics in DGP 
obeys the usual 4D diffeomorphism invariance, we can al- 
ways pick a very fast-moving observer to explore the field 
of a massive source, and transform all of the relevant 
physics to her rest frame. In her frame, thanks to Rela- 
tivity Principle, the source mass would appear to move 
with a very high speed v, and therefore its gravitational 
field should be well approximated by our shock wave so- 
lutions. For a sufficiently fast observer, the corrections 
due to th e mass s hould be suppressed by the powers of 
m/p — y^l/v 2 — 1. This suggests that if organized as an 
expansion in powers of m/p, perturbation theory may be 
under control (however, by itself this may not help with 
any ghost infestations on the self-inflating branch [3). 

The paper is organized as follows. In the next sec- 
tion we construct the background solutions for a ten- 
sional brane without particle excitations on it. They are 
de Sitter vacua of 5D DGP. In section III we present a 
detailed derivation of the gravitational shock waves, and 
show that they can be written in terms of two-variable 
hypergeometric functions. We outline a spectacular new 
non-perturbative channel for production of bulk gravi- 
tons on the self-inflating branch, when gravity is modi- 
fied at exactly the brane de Sitter radius, in section IV, 
and consider the limiting form of the solutions at short 
distances, as well as several special cases in section V. 



We then turn to the problem of gravitational scattering 
of relativistic particles in section VI, and summarize in 
section VII. 



II. DE SITTER VACUA 

Brane-induced gravity models Q are given by a bulk 
action with metric kinetic terms in both the bulk and on 
the brane. In the case of one extra dimension, 

S = J d 5 x^E^-R 5 - J tfxy/giMlK 

// M 2 \ 
^Xy/gi(-^- i? 4 - A - £ mat tcr) • (1) 

Here gAB and g M „ = d^X A d v X B qab are the 5D metric 
and the induced 4D metric on the brane, respectively, 
i?5 and i?4 are their Ricci tensors, K — g AB Kab is the 
usual Gibbons-Hawking term necessary in space-times 
with boundaries (defined as the trace of the extrinsic 
curvature of the brane, for variational principles see e.g. 
pCiL and A and £ ma tter the brane tension and mat- 
ter Lagrangian (which we separate from each other ex- 
plicitly). The indices {A,B} and {/Lt, v\ are bulk and 
3-brane world-volume indices, and the gravitational cou- 
plings are controlled by the bulk and brane Planck scales, 
M4 and M5, respectively. 

Varying this action away from the brane one finds the 
gravitational field equations in the bulk, which in the case 
of empty bulk admit locally 5D Minkowski solution. The 
brane dynamics can then be incorporated via the modi- 
fied Israel junction condition [Til l2fij. which in the case 
of brane-induced gravity in addition to the brane extrin- 
sic curvature also receives corrections from the intrinsic 
curvature as well, arising from the variation of -R4. This 
junction condition can be viewed as the equation of mo- 
tion of the brane, defining the world- volume the brane 
sweeps in the locally Minkowski bulk. In the case of 
empty bulk, one can readily transform to the rest frame 
of the brane by going to the Gaussian normal coordinates 
for the 5D metric, ds\ = dw 2 + g^{w, x)dx tl dx u , and in- 
corporate the brane equations of motion as (^-function 
terms in the bulk field equations, taking a very compact 
form. In 5D, those field equations are p| 

MlGi B + MiG^5 A S B 5(w) = -T\5 A 6 B 5{w) , (2) 

where G A b, G^ v are bulk and 3-brane-localized Einstein 
tensors, and T^ v is the brane stress-energy, respectively. 
In this coordinate system one can easily orbifold the bulk, 
by imposing a Z2 symmetry w — > — w around the brane 
situated at w — 0. 

For a brane with tension A 7^ and without any local- 
ized or distributed matter sources, T M „ = ~\8^ v . The 
corresponding solutions describe de Sitter vacua on the 
brane. The symmetries of the brane-bulk system then 
dictate the form of the metric, which must be given by 
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the warped 4D de Sitter space in a locally flat bulk |25j 
dsl = {l-eH\w\fdsl dS + dw 2 . (3) 



The replacement w — > \w\ implements the orbifolding 
procedure. Aside from the parameter e, which can take 
values ±1 in this metric is identical to the 5D version 
|2m of the inflating Vilenkin-Ipscr-Sikivic domain wall 
[29| . In what follows we shall use the static patch de 
Sitter metric, 



ds\ds 



(1 - H 2 r 2 )dt 



dr 2 



(1 - Hh 



r 2 dfl 2 . 



(4) 



instead of the spatially flat one utilized in 28]. The 
space-time of © is a 4D de Sitter hyperboloid in a 5D 
Minkowski bulk. Here e arises because we can retain ei- 
ther the interior e = +1 or the exterior e = — 1 of the 
hyperboloid after orbifolding, thanks to the brane intrin- 
sic curvature R4. The junction conditions relate the 4D 
curvature and the brane tension as 1251 



H 



2M( 
' Mi 



H 



A 



3M| 



(5) 



The simplest way to derive this equation is to substitute 
the metric © into the field equations @, recall that 
away from the brane is just a different parameteriza- 
tion of the locally flat Minkowski metric, and compute 
the terms proportional to S(w), which in this case arise 
from the double iu-derivatives in the bulk terms on the 
LHS in J2J, and from the brane curvature terms. Then 
matching them with the tension source yields J3J). Wc 
will revisit this derivation in more detail in section III. 

Focusing on the solutions with H > (because H < 
cases are their PT transforms), on the normal branch 
defined by the choice e = +1, we see that the equation 
|J5J reduces to the 4D Friedman equation in the limit 
M 5 < M 4 , 3H 2 ~ A/Mf. In this case in the bulk we 
keep the interior of the hyperboloid. It has finite volume 
and so the perturbative 4D graviton exists. The brane 
curvature terms ~ Mj J d A x^fglRn/2 suppress the cou- 
plings of the m g > graviton KK modes, producing the 
4D effective theory when M5 <C M4. In Fig. ^ this case 
corresponds to keeping the unshaded region of the bulk, 
i.e. the interior of the de Sitter hyperboloid. 

On the other hand, on the self-inflating branch, defined 
by e = —I, at low tensions A <C I2Aff /M|, from the eq. 
© one finds H ~ 2Mf/Mf. In this case in the bulk 
one keeps the exterior of the hyperboloid that opens up 
to infinity, corresponding to the shaded region in Fig. ^ 
This bulk volume is infinite, and so there is no perturba- 
tive 4D graviton. The effective theory arises only from 
the exchange of the bulk resonance. 

We note that this solution resembles Witten's "bubble 
of nothing" [3(j (see also [U, see Fig. [I]), which first 
collapses from infinite size to a minimum, and bounces 
back to infinity, gobbling up the bulk along the way. 
However now the bubble wall carries finite energy den- 
sity. This may help protect the bulk from being "eaten" 



Static 
patch 



de Sitter 
brane 




FIG. 1: Penrose diagram of a de Sitter brane in a flat 5D bulk. 
We suppress one angular coordinate on S 2 on the brane and 
depict another, so that the diagram looks like two cones glued 
along their base, with the hyperboloid of the brane world 
volume inscribed inside the cones. In the case of the normal 
branch, in the bulk one retains the interior of the brane world- 
volume, denoted with e = +1. For the self-inflating branch, 
one retains the exterior, marked e = — 1. The static patch on 
the brane is the interior of the diamond-shaped region on the 
brane world- volume. 



by "nothing" too quickly. Namely, the probability of 
nucleating a "bubble of nothing" is given by the expo- 
nential of the negative action of its Euclidean extension, 
V ~ exp(— Se)- If we take the self-inflating branch so- 
lutions (j3J and glue them to a Euclidean hemisphere at 
the minimum radius to model bubble nucleation, we find 
that the corresponding action is 



^(A + 16Mfff) 



(6) 



where V4 is the volume of the sphere of radius H 1 (for 
a simple calculation of the action with Gibbons-Hawking 
terms see e.g. [H Ej^). In the limit A < I2Mf/Mf , 
M s <C M4, only big bubbles with a large Euclidean ac- 
tion, Se ~ (M4/M5) 6 3> 1, are nucleated. Therefore, if 
the tension and the bulk Planck scale are somehow fixed 
to be much smaller than the 4D Planck scale M4 (see e.g. 
[33j for the discussion of scales in brane-induced gravity) 
the nucleation rate of "bubbles of nothing" may be suffi- 
ciently slow to render the self-inflating branch long-lived. 



III. DERIVATION OF SHOCK WAVES 

We now turn to the construction of the gravitational 
shock wave fields generated byrelativistic particles mov- 
ing on the brane in DGP [15|. Imagine that a pho- 
ton is moving on a tensional brane, with a momentum 
p = 27TJA Because it is a relativistic particle trapped on 
the brane, it moves along a null geodesic of the brane- 
induced metric. By the Principle of Equivalence, its mo- 
mentum generates gravitational field, since it contributes 
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to the stress-energy tensor. There are two approaches for 
generating the gravitational field solutions. One could 
start with a linearized solution for a massive particle at 
rest, and boost it to relativistic speed, simultaneously 
taking the limit m — > such that m cosh 7 = p remains 
constant. This is how the solutions in 4D were originally 
found 0, ^3 ■ The procedure is completely analogous to 
the method of generating the electro-magnetic potential 
of a charged massless particle (for a nice review, see pi)). 
The same procedure has also been applied to the massless 
particles moving in maximally symmetric spaces |22t I23I ] 
and Randall-Sundrum braneworld models |3J,|3a|. The 
reason why this procedure works as a method of generat- 
ing exact solutions of gravitational field equations is that 
if we expand the exact metric for a massive source in pow- 
ers of the mass m, after a boost to relativistic speeds the 
quadratic and higher powers of m vanish because there 
is only one factor of the boost parameter cosh 7 in the 
metric, which is compensated by one power of m when 
COSI17 — * 00, m — > and mcosh7 = p = const. Hence 
the linearized solution becomes exact. Essentially, in the 
extreme relativistic limit the solution becomes scale in- 
variant, and the non-linear effects which are controlled 
by the fixed ratio 171/M4 drop out. 

The previous discussion on boosting linearized solu- 
tions describing point masses to generate exact solu- 
tions of gravitational field equations illuminates the na- 
ture of the solution as a planar shock wave. Because 
of the infinite boosting the field lines of the gravita- 
tional field become confined to a plane perpendicular 
to the direction of motion. In terms of the null coor- 
dinate u which parameterizes the world-line of the par- 
ticle v, the field experiences a jump at u = encod- 
ing causal propagation of the source. When the parti- 
cle flies by an observer, its gravitational field perturbs 
the observer at precisely that instant 0. This prop- 
erty allowed Dray and 't Hooft to develop a very elegant 
cut-and-paste technique [18| to find solutions, which they 
applied to asymptotically flat 4D backgrounds. Higher- 
dimensional examples were considered in IHI23. This 
technique was later applied to general 4D GR back- 
grounds by Sfetsos [24| . To circumvent the need for a de- 
tailed form of the linearized solution for a localized mass 
in DGP, we will apply a variant of this cut-and-paste 
technique to DGP. We start by rewriting © in suit- 
able null coordinates, defining u = 77 \/ 1+^1- exp(Ht) 



and v = -jj 



exp(— Hi). We also introduce a new 

bulk coordinate \z\ = — ^jTn(l — eH\w\). It simplifies 
evaluating the curvature tensors with a conformal map 
as in [2a|. Note that S(w) = S(z) because the warp fac- 
tor and its derivative are equal to unity on the brane. In 
terms of these coordinates, we can rewrite the metric Q 
in the following form: 



eH\z\(_ 
1(1 



Adudv 



(1 + H 2 uv)' 



■ + ( 



H uv 2 dfl 2 



1 + H 2 i 



H 2 



(7) 



Suppose now that a photon moves along a null geodesic 
on the brane, let's say the ii-axis of the coordinates uti- 
lized in (0. The equation of this null geodesic is u = 0. 
If we choose the origin of this photon at infinite past as 
the North Pole of the static patch, an observer located 
there sees a photon streaming away along the past hori- 
zon. Following Dray and 't Hooft, we introduce the shock 
wave on top of the background JJJ by including a jump 
in the v coordinate at u = 1181 . We substitute 



v — 
dv 



v + Q(u)f, 
-> dv + 0{u)df, 



(8) 



where / is the shock wave profile, whose precise form is 
to be determined by requiring that combining J7J and 
||SJ| still remains a solution of the field equations. The 
wave profile / depends only on the spatial coordinates 
transverse to the wave, which in our case are the angles on 
the 2-sphere and the bulk coordinate z. Here 6(w) is the 
Heaviside step function. Note that in flat 4D background 
© is a trivial diffeomorphism if / = const. [18|. In that 
surroundings, if df 7^ 0, the replacement rule for dv in © 
ensures this is not just a diffeomorphism. Changing the 
coordinates to v = v + Q(u)f we find dv — > dv — 5(u)fdu. 
Substituting v, dv — > v,dv — S(u)fdu in J7J and dropping 
the carets, we obtain the ansatz for the metric with the 
shock wave included: 



dsj 



-2eH\z\( 

1(1 



Adudv 



(1 + H 2 uv) 2 

,1-H 2 uv, 
+ ( 



4S(u)fdu 2 
(l + H 2 uv) 2 
dflo 



l + H 2 i 



Note that with this, the ansatz for the induced metric on 
the brane is 



ds\ 



Adudv 



A5{u)fdu 2 
(l + H 2 uv) 2 ~ (l + H 2 uv)'- 
,1 



+ 



H 2 uv 2 dfl 2 
1 + H 2 uv> ll 2 



(10) 



Now we take the ansatz (JSJ, (|1U|) and substitute it into 
the field equations (J5J, where we also include the pho- 
ton stress-energy contribution in the source on the RHS. 
Thus the total stress-energy tensor with a photon moving 
along u = is 

T\ = -XS" U + 2^=g4 UV 5(e)6(<p)d(u)S%6» . (11) 



We have chosen the coordinates on the 2-sphere such that 
the photon trajectory is at 9 = 0. 

To evaluate the curvature we utilize the conformal 
transformation trick as in |3(j, splitting the metric as 
ds\ = fl 2 dsl, where Q = exp(— eH \z\), and 



dsi 



Adudv 



4S(u)fdu 2 

(l + H 2 uv) 2 ~ (1 + H 2 uv) : 



+ 



, 1 — H 2 uv n9 dQ,o , o 
1 + H z uv H z 



(12) 
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A straightforward albeit tedious computation, where we 
treat S(u) and its derivatives as distributions, and use 
the distributional identities uS(u) = 0, u 2 S 2 (u) = and 
f(u)5'(u) = -f'(u)6(u) EE!, yields the Ricci tensor 
for the conformal metric (|12l) . The stress-energy is 
automatically conserved, V^T^ oc u8(u) = 0, so the 
matter sector field equations are satisfied. We note that 
this computation also yields automatically the compo- 
nents of the Ricci tensor of the induced metric on the 
brane, which appear as the subset of the curvature ten- 
sor for the conformal metric (|12[) . because f2(0) = 1 and 
the coordinates in are Gaussian, as is obvious from 
comparing (|10|l and (|12() . The resulting non- vanishing 
Ricci tensor components are 

R 5 uu = 25(u)(d 2 f + H 2 (A 2 f~4f)), 

R iuu = 2H 2 S(u)(A 2 f-4f), 

R§ uv — ^4 UV 



6H 2 
(1 + H 2 i 



R$ab — Riab — 3 



(t 



H 2 UV\2 



H 2 i 



9ab , 



(13) 



where R§ ab and R4 M „ are the Ricci tensors of the metrics 
(I12J1 and l)10[l. and a, b and g a b are the indices and the 
metric on the unit 2-sphere, respectively. The operator 
A2 is the Laplacian on the unit 2-sphere. 

Using the conformal transformation ds\ — £l 2 ds\ 
which gives the relationship 

R 5 ab = Rsab -3V^V B lnn + 3VAlnOV B 

- 9AB (V 2 In fl + 3(V In O) 2 ) , (14) 

we finally find 

R 5zz = 8eHS(z), 

R 5uu = 25(u)(d 2 J - ZeHd [z[ f + H 2 (A 2 f + 2/)) 

-8eHfS(u)S{z) , 
AeH 



R 



5 uv 



Rhab — T7 



(1 + H 2 uv) 2 
2e / l - H 2 uv\i 



S(z)g a b 



(15) 



H\l + H 2 uv; 
Although Ricci tensor vanishes away from the brane at 
z = even when the shock wave is present (-R5 uu in H15|) 
vanishes at z ^ by virtue of the field equations, see 
below), the bulk is not Minkowski any more, since the 
waves extend off the brane and deform the bulk. One 
would find their non-vanishing contributions away from 
the brane in the Riemann tensor. Now from (|13fl and l|15|) 
we can obtain the explicit form of the field equations J2J- 
First, it is easy to check that we can rewrite @ as 



MlR 5 



M 2 4 
-fRi5{z) = -\8(z) 

1 



(16) 



M£R 5llv + Mi{R ijxv - -gi^Ri)5{z) 



Here S(Q) is the covariant d- function on the unit 2-sphere, 
6(fl) = 5{9)5{<j>) I 'Jgi = 5(cos9- l)S(<j>), which peaks at 
9 = 0. Substituting O and (JTHJl in the first of itTHI) 
we recover just the background equation JSJ), linking the 
Hubble parameter H and the brane tension A. The new 
equation for the wave profile / comes from the uu com- 
ponent of H16(l . Indeed, direct substitution of (|13(l and 
(|15fl yields a linear field equation for the wave profile: 



M 2 H 2 



(d 2 J-3eHd lzl f + H 2 (A 2 f + 2f) 



+ (A 2 f + 2f)8(z) = ^S(Q)S(z).(17) 

The remaining components of the second of (|16fl again 
reduce to the background equation (JSJ. Thus ((T7Jl fully 
determines the shock wave profile /. Note that in the 
limit M5 — > the bulk derivatives disappear, and we 
can factor out S(z) to recover the correct 4D de Sitter 
equation of jU 0] • 

Now we can turn to solving H17|) . It turns out that one 
has to carefully implement the boundary conditions for 
the single particle problem (|17fl in order to avoid addi- 
tional singularities that would arise because the 2-sphere 
on which the particle propagates is compact [24[. These 
boundary conditions are handled in the simplest way by 
adding an extra particle with the same momentum p, 
run ning in the opposite direction in the static patch Q 
[22I l23j | In fact, this two-source solution correctly rep- 
resents the limit of infinite boost of the Schwarzschild-dc 
Sitter geometry plE^. So following [illHIH we add 
an extra term, -$z8(Vl')8(z), on the RHS of JTiJ. Here 

5(0') is peaked on the opposite pole on the 2-sphere, at 
9 = 7r (see Fig. [2J. Then we can return to a single source 
by taking this solution and multiplying it by Q(it/2 
as in 




9 = 



rS4^(z) + H 2 p(g iuv ) 2 5{n)5{u)5{z)5^5 u v 



North Pole, r = 

FIG. 2: Trajectories of two photons in the reference frame 
of an observer at rest on the North Pole, at r = in the 
static patch. The photons move along opposite directions on 
the S 2 , at coordinate distance ru = l/H from the pole. We 
are explicitly depicting the polar angle 8 and suppressing the 
azimuthal angle <j> on S 2 . 
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So let us consider eq. I|17|l with two photons, when one 
of them moves along 9 — and the other along 9 = tt: 



Mf 
M?H 2 



(d 2 J-3eHd lzl f + H 2 (A 2 f + 2f) 



+ (A 2 / + 2f)6(z) = Jfc (<5(fi) + 5(0')) S(z) . (18) 

Since the sources on the RHS are S(£l) = S(cos9 — l)6((f)) 
and 5(0') = 5(cos9 + l)5(<fi), they can be decomposed in 
terms of spherical harmonics as 



z=o m=—l 

oo Z 

j(n') = E E ^(0,0)^-0^). (19) 

;=0 m=-i 

Since we chose the photon trajectories to be along 9 — 
and 9 = 7r, we can use the addition theorem for spherical 
harmonics to replace them with Legendre polynomials 
P;(cos£?) in the expansion for the wave profile. Indeed, 

from EL-i^fO.O)^^,^) = ^iPi(cosfl) and the 
similar equation for the 9 — > tt — 9 terms, and linearity 
of 118|) , we deduce that the solutions are of the form 

oo 

/ = E(// +) w^( cos0 ) + ^ ) (^(- cos0 )) • ( 20 ) 

1=0 

Here // ±} (z) are the bulk wave functions; /, is sourced 

by the photon at 6 = and /, ' by the photon at 9 = tt. 
By orthogonality and completeness of Legendre poly- 
nomials, using the expansions for the 5- functions from 
above, the field equation 1|18[) yields the same differential 
equation for the modes ft (z): 

d 1 J l - teHdufi + H 2 (2 - 1(1 + I)) ft = 



M 2 H 2 ( (21 + l)p 



Ml 



2ttMI 



(2-l(l + l))f l )6(z). (21) 



We should interpret 5-function on the RHS by pillbox in- 
tegration as a matching condition for the first derivatives 
of on the brane, recalling also the orbifold symmetry 
which enforces fi(—z) = fi(z) . This yields the boundary 
value problem for // on the brane: 



d 2 Jt - 3eHd z ft + H 2 (2 - 1(1 + 1))/, = , 
fl(-z) = f t (z) , 

2-1(1 + 1) H21 + 1 



(22) 



-Hh(0) = 



g 2ttM( 



P21 (cos 9) , circumventing the unphysical 4D singularities 
in / = 1 terms like in pll^: 



f = 2j2f2i(z)P2i(cos9) 



(23) 



1=0 



We still need to specify the boundary conditions for / 
in the limit z — > 00. We do so by requiring square inte- 
grability in the bulk, so that the solutions are localized 
on the brane. To see how to implement square integrabil- 
ity on the modes /; , it is sufficient to consider the simpler 
example of a free massless bulk scalar field, with action 

S = - J d 5 x^(\7 x ) 2 . (24) 

In terms of the conformal geometry (|12J) . the action is 

S = - J dzdtxy/gi n 3 (Vx) 2 , (25) 

where O = exp(—eH\z\). Therefore, the norm for the 
expansion of x i n terms of the 4D modes satisfies 



| Xq || 2 cx / dze-^ X 2 a 



(26) 



We are interested in the modes of \ which are simple 
exponentials of \z\, because from (|22|l it follows that /;'s 
are such simple exponentials. Hence let Xa(x) oc e~"' 2 '; 
then from H26|) we conclude that a mode will be nor- 
malizable only if a > (1 — 3eH)/2. Solving the differ- 
ential equation in l|22(l we see that the modes are of the 
form ~ e [±2i-(Ti-3e)/2] H\z\ _ jj ence we should retain only 
the sub-leading bulk mode which is square-integrable, or 
equivalently, localized on the brane: 



ft = A,e" [2 ' +(1 ~ 3e)/2IH|z| 



(27) 



We can now write down explicitly the series solution 
of (|22|l which satisfies these boundary conditions. Sub- 
stituting (|2*7j) in the last of (|2*2"|) and solving the ensuing 
algebraic equation for we find 



/(«,*) = - 



2ttM- a 



E 



4Z + 1 



1(21-1 + 



)(l + l + 



(H 



1=0 V" 1 2 A* 1 * 4 

xe -[2l+(l-3 £ )/2]»l»lp 2I ( cosfl ) j (28) 



As a check, in the limit g = 0, at z = this reduces 
to the 4D series solution of for both e = ±1. We 
can write the series (|28|) in the integral form. Defining 
r = e~ R \ z \ and x — cos 9 and factorizing the coefficients 
of the expansion, we first find 

[3-(l-e)g] P ^ T 2i + (l-3 e )/2 



where g = 2Ml/(M 2 H) = l/(Hr c ) (see Q). Since both *) = ~ (3 + es]nM 2 E 



/, solve the same boundary value problem, /, - 

ft = fi- Because Pi(— x) = (—l) l Pi(x), the solution 
will be an expansion only in even-indexed polynomials 



-P2l(x) 



[3+(l + e) 8 ]p 
2(3 + eg)TTM 2 j^l + 1 



E 



r 2i + (l-3e)/2 



(l+<0l 



-P 2l (x) . (29) 
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In the limit g = 1 the I = term in the first series 
diverges for e = — 1. There are no divergences in the 
case g = 3, as we see from eq. We postpone the 

discussion of the limit g = 1 until the next section, and 
for the remainder of this section take g ^ 1. 

Recalling the definition of the generating function for 
Legendre polynomials, (1— 2xt+t 2 )~ 1 / 2 — ^2^q Pi(x)t 1 , 
and using 



lim / dd$ l+0 



r 



Z+/3+1 



„ lim 

1 + + 1 e^O 



dM l+p (30) 



to regulate divergences when I + (3 < —1, we integrate 
the generating function expansion over r. This yields 
the identity 



™ T 2l+0+l 



P2l{x) 



(31) 



- lim (2 

2 e^oI 



«<-(l+/8)/2 



1=0 



VI - 2xd + d 2 VI + IxV + i? 2 



Applying eq. I|31|) to eq. i|29|) , after straightforward ma- 
nipulations and performing explicitly the limit c^Owe 
find 



m,z) 



p 



l-3e 
T 2 



P 



2ttM 2 1 - ±^g 2(3 + eg)7rM| 



1 



1 



/ dd( , = 

lo VVl-2xtf + d 2 Vl + 2x7? + d 2 



(32) 



[3-(l-e)gk9 



5(1-0-4 



[3 + (l + e )g]^ 

(g+3)(l + Q 



When g V 1 the integral (|32ll is finite and well-defined 
everywhere except at x = ±1, where it has the usual 
short-distance singularities, because it is a Green's func- 
tion in the transverse spatial directions and x = ±1 is 
where the sources are located. We remove the spurious 
poles in the limit e — > in the integral representation 
of (|32|) , arising from the implementation of the regulator 
(|30|l . I|31l) . by the introduction of the factor —2 in the 
second line of the integrand of (|32l) . 

For completeness' sake, we note that the integral (|32f) 
can be expressed in terms of the two-variable hypergeo- 
metric functions Fi(a,/3i,/32,j;p,q) |3JJ, which are de- 
fined by the integral representation 



F 1 (a,0i,02,a + n;p, q) 



d( x 



x r-^i - o M_1 (i -Kr^a - qcr 02 m 

valid when Re a > 0, Re /1 > 0. As an illustration, we will 
only consider here the shock wave profile / on the brane, 
at z — (i.e. r = 1), on the self-inflating branch e = —1. 
Similar derivation can also be performed for the normal 



branch e = 1. We will split up the integral in 132(1 into a 
sum of two integrals, with terms oc (l + 2x , d + 'd 2 )~ 1 / 2 — 1, 
and evaluate each integral separately by using the Eu- 
ler substitutions 1 + ^2(1 + x)( = Vl T 2xtf + d 2 + -3. 
To see how the hypergeometric functions in l|33|) emerge, 
consider e.g. 



dti( , 

Vl - 2xd + t) 2 



1)(V^^ 2 + T^^' 
3-g 3-g 

(34) 

and plug in 1 + ^2(1 - x)( = Vl - 2xd + d 2 + After 

straightforward algebra, we can rewrite the integral 13411 , 

after introducing Q = —y/(l — x)/2 = — sin(0/2), as 



3-g 



g-2 



g ^ C-Q 
3-g C-Q - 1 

Breaking up the integrand into the powers of the kind 
that enter in 1|33|) we can rewrite (|35|l as 



I = 



(3 - 2g)r(g - 2) xFx(g - 2, 2 - g, g, g - 1; Q, i) 



(3-g)r(g-l) (-Q)s 

(3 - 2g)r(g - 1) Ms - 1, 2 - g, g, g; Q, £ 



+ 



(3-g)r( g ) (-g)s-i 

3 2^(1, -1, 3, 2; Q,^) 



(3-g) (-Q) 3 

3 iq(2,-l,3,3;Q,i) 



(36) 



2(3 -g) (-Q) 2 

The remaining integrals in (|32|l can be found from i|36|l by 
substituting x — ► —x. Then, using P = — <J (1 + x)/2 = 
— cos(0/2), we finally obtain 



m) = 



2irM 2 1 - g 2^(3 - gi.\/ : 



(3-2g)r(g-2) sFi(g-2,2-g,g,g- 1;Q,^) 

r(g-i) (-Q)« 

(3 - 2g)r(g - 1) Fi(g - 1, 2 - g, g, g; Q, |) 

r(g) {-QY- 1 

1*1(1,-1, 3,2; Q,£) Fi(2,-l,3,3;Q,i; 



(-Q) 3 2(-Q) 2 
(3 - 2g)r(g - 2) xFtjg - 2, 2 - g, g, g - 1; P, ^ 



r(g-i) (-p)b 

(3-2g)r(g-l) Fife-M-g.g.gjP, £) 



r(g) (-P)s- 1 

a;Fi(l,-l,3,2;P,^) Fi(2, -1, 3, 3; P, ^ 



(-P)E 



2(-P)= 



(37) 

')}■ 



In the expressions for the hypergeometric functions one 
would encounter an additional, spurious, logarithmic sin- 
gularity. This singularity is manifest in the integrals l|34|) 
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or l|35|) in the lower limit of integration. However, after 
Euler substitutions for integration variables, one has to 
carefully relate the singularities in (|3"4*j) and fl3"K)l. recall- 
ing that these integrals are defined as principal values of 
integration over in accordance with the regularization 
prescription in H3()|) . After this is done, these spurious 
logarithmic singularities cancel precisely between the Q 
and P dependent terms in (|37|) , leaving only the physical 
short distance singularities carried by the null sources on 
the brane. We will revisit this more closely in section V. 



IV. RESONANCE ON THE SELF-INFLATING 
BRANCH 

As we have mentioned above, the solution (|29|l for the 
shock wave displays a spectacular behavior on the self- 
inflating branch, e = — 1, asg— > 1, In this limit the 
tension vanishes, as can be seen from eq. (|SJ). The I = 
term of the first sum has a pole there. No such poles arise 
on the self-inflating branch, where the series is finite. 
What then does this divergence mean? To gain some 
insight we revisit the boundary value problem 122|) for 
the I = mode but slightly relax the boundary conditions 
in the bulk, by retaining the / = solution which is not 
square-integrable as well. Writing this mode as 



fo = A e- 2H ^ + B e 



-HI. 



(38) 



and solving for Aq and Bq we find that when g = 1 the 
mode oc A is an eigenmode of the homogeneous (p = 0) 
boundary value problem l|22l) . completely cancelling on 
the LHS of the last of eq. (|22l) . This is the origin of 
the divergence in l|29|l for g = 1 and mathematically is 
an example of the Fredholm's alternative: either the in- 
homogeneous differential equation or the associated ho- 
mogeneous differential equation (i.e. with non-vanishing 
source or with vanishing source, respectively) have a so- 
lution, but not both at the same time. On the other 
hand, the mode oc Bq remains. Hence we could solve the 
last of (EU for it: 



Bq = 



2nM 2 



(39) 



Because we have excluded the non-integrable bulk modes 
with I > 0, this mode is the fastest-growing one in the 
limit z — ► oo among all the allowed solutions of (|22|l for 
any I. Thus far from the brane the resulting wave profile 
behaves as 



V 



TtM 2 



-H\z\ 



(40) 



Hence as z approaches the Cauchy horizon at z 
(see Fig. ^) the asymptotic metric © is 



eH\z\( 

1(1 



Adudv 1 — H 2 uv, 2 dfl2 

(1+H 2 uv) 2 + A + H 2 uv' IF 



dz 2 



Ap 



e-' H ^6{u)du 2 . (41) 



Because the last term remains finite, the shock wave is 
not localized to the brane, but extends deep into the 
bulk, as anticipated. This solution carries a lot of en- 
ergy and curvature into the bulk, similarly to the grav- 
itational pp-waves in the Randall- Sundrum setup, stud- 
ied in 38]. Indeed, we can rewrite the asymptotic metric 
(|41|l expressed in coordinates u, v, z in terms of the orig- 
inal static patch coordinates r, t, w (see the text above 
eq. J7J) ) , and transform further to the spherical polar 
coordinates in the bulk far from the brane, defined by 

1 

H 

1 P 



: sin if , 



sinh(iJi) = 



\J p 2 cos 2 tfi 



(42) 



Then the asymptotic metric (|41f) becomes 



ds 2 5 



-dr 2 +dp 2 +p 2 dn 3 



8pVP 2 



■kMIH 



-5(X 2 -l)dX 2 



where X = 



p COS ip + T 



(43) 

— , and dfls — dip 2 + sin 2 (pdQ,2 is 

v/p 2 — t 2 -fp sin ip 

the line element on a unit 3-sphere. Since this asymptotic 
metric is invariant under S 2 rotations, the solution looks 
like an imploding shell of energy in Minkowski space. For 
a hint of how such a shell is generated, move away from 
g = 1 and reconsider the boundary value problem i|22|) 
while still keeping the non-integrable bulk wave function 
in the / = mode, as in eq. (|38|l . Then 122|) relates Aq 
and Bq according to 

2^(^k- 2{i - £}A °y (44) 

Thus we see that when g ^ 1, the localized mode oc Aq 
can absorb fully the effect of the source oc p. This is 
energetically favorable, because the non-integrable mode 
oc Bq decouples because it has vanishing overlap with 
the brane where the source is located. Hence we can 
set Bq = and forget about it. However in the limit 
g = 1 the modifications of gravity occur at the scale equal 
to the cosmological horizon on the brane. Brane and 
bulk begin to "resonate" : in a time-dependent problem, 
as g — > 1 slowly, the localized I = mode would grow 
ever larger, spreading further out into the bulk. It would 
develop a larger overlap with the bulk graviton, so that 
it can produce delocalized gravitons in order to shed the 
energy which drives its unbounded growth. In this way as 
g — * 1 energy could escape into the bulk. This indicates 
an onset of a dramatic new instability, which warrants 
further investigation. 



V. SHOCK WAVES AT SHORT DISTANCES 

In this section we consider / on the brane, when z = 
(i.e. t — 1), at transverse distances TZ well inside the cos- 
mological horizon, TZ <§; H" 1 . Note that we can rewrite 
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eq. 10) as (l + eg)H 2 = A/(3Mf). This shows that on the 
normal branch g can vary between zero, in which case we 
recover the 4D limit, and infinity, when we recover the 
5D limit, with the "resonance" discussed in detail in the 
previous section at g = 1. On the other hand, on the self- 
inflating branch, g < 1 or r c > 1/H, as long as A > 0. In 
the limit g — > 1 on the self-inflating brane, as mentioned 
before, the tension vanishes. 

Let us start with g = 0. In this case ll-i'2t is identical to 
the 4D case on both branches. Indeed, it can be expressed 



as / 



p 

27ri\/r 



(1 — I_ — J_|_) where are the integrals 



VI T 2xd + $ 2 



1 



(45) 



which obey I+(x) = I—(—x). Hence as in section III 
it is enough to evaluate one, say J_. Using the Eulcr 
substitution £ = Vl ~ Zxtf + tf 2 + d we can rewrite the 
integral 7_ as 



I- = 4x 



l+s- 



(£+l) 2 (£-l) 



(e + i) 2 



1 f 1+ V^^ (z - l) 2 (z + 1) 

1+x [i+V^ (z-l) 2 (z + l) 

where in the first line we have replaced 1 — > 1 + £_ in 
the lower limit of integration in order to regulate a spu- 
rious logarithmic singularity, which we have discussed in 
section III. A direct evaluation of the integrals yields 



T 1 2 1 

1- = x m x 

s- 2 



(47) 



and thus, using I+(x) — x), and noting that the 
regulator is now e+, also 



, 2 1 

7+ = —x In , 

e+ 2 

such that the final result for I = I + + i_ is 



I = x In — - 1 . 

£_ 



(48) 



(49) 



To determine the ratio of the regulators recall that 
they are the image of the regulator e for the variable ■& 
under the Euler substitution for the integration variables 
in I±. That yields e T — (1 =p x ) £ , so that 



e+ 1 + a: 
e_ 1 — x ' 

Therefore the singularities indeed cancel out, and 

/ = x l ri [l±^l_l. 
1 — x 



(50) 



(51) 



Substituting this in the equation for the shock wave pro- 
file finally yields 



P 



2ttM, 



2 Mr 



(52) 



The metric transverse to the null particle on the brane 
is g?s 2 U=m=o = (HI2/H 2 . Therefore the proper "radial" 
transverse distance for small polar angles 9 is 1Z ~ #/-ff , 
so that a; = 1 — H 2 1Z 2 /2. Thus the shock wave profile 
(|52|l reduces precisely to the flat 4D solution [l7l llq : up 
to OJJZ 2 1 H~ 2 ) corrections, using the sign conventions of 
Ell, we find 



p 



71 



7rM| 



ttM 2 



(53) 



As long as it is finite, the numerical value of the constant 
term in this equation is physically unimportant since a 
finite constant term in the shock wave profile in a flat 
4D spacetime can always be changed by a simple 4D dif- 
feomorphism [l^. We cannot write down the integrals 
which incorporate the 0(g) corrections in closed form 
because they contain terms like J G?£ln(l + &0/C H3- 
Nevertheless, it is easy to extract their short-distance ef- 
fects: at short distances the singular 0(g) powers in the 
integrand precisely cancel out. In the remainder, since 
the leading order singularity in l-il'li is logarithmic around 
g = 0, the additional logarithms soften the corrections 
further, and render them finite as 1Z ~ > 0. Hence at 
short distances the solution looks exactly the same as in 
4D, and the corrections are negligible when 1Z < H^ 1 /g. 

In fact this behavior persists for any finite value of g. 
We can see this from identifying the leading contribution 
in (|32f> at short distances. Rewriting the integrals (|32f> by 
using the same substitutions £ = Vl T 2x$ + 1) 2 + $ and 
looking at the limit 1 — x <C 1, we identify the leading 
singularity. After cancelling the spurious terms as above, 
for all values of g 7^ 1 and e the leading singularity at 
short distances comes from the term 



P 



l+V2(X-a:) 



P 



nM 2 



In 



1 



+ finite . 
( 54 ) 

Substituting x = 1 — H 2 lZ 2 /2, we recover the leading 
logarithm in (|53|l . The corrections to this result at larger 
distances in general differ from the corrections which 
arise when one expands the 4D de Sitter solution 152|) 
around the short distance limit describing shock waves 
in asymptotically flat space (|53|l . where one finds correc- 
tions to come as even powers of H1Z. In DGP, however, 
the corrections in (|55|) will arise as odd powers of H1Z 
as well, starting with the linear term. This signals the 
presence of the fifth dimension concealed by the "gravity 
filter" . Therefore in general we will have 



((l+ aij ff2 7e2 + ... )ln( __ ) 

+const + ^HTZ + b 2 H 2 K 2 + ..) , (55) 



irM 2 
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where the coefficients , bk are numbers that can be com- 
puted explicitly for given values of g and e, and we have 
seen that e.g. b\ ~ g. 

To see how 5D gravity reemerges, we can take the limit 
g — ► oo on the normal branch e = 1, keeping M5 fixed. 
This removes the brane-localized terms oc M4 from the 
action and neutralizes the "gravity filter" . In this case 
the background solution reduces to the inflating brane 
in 5D Minkowski bulk |28||. From eq. JUJ) we see that 
the shock wave profile on the brane at z — in the limit 
g — ► 00 reduces to 



/ 5D (0) = 



ttM 5 3 ^ 21 K ' 2ttM| ^ 21 - 1 ' 



(56) 



2=0 2=0 

where we have used Mfg = 2M§/H. We can evaluate 
the first sum using the identity 



TP2i(x) = 



1 



1 



2=0 



(57) 



V 8 (! - *) V8(l + a:) 
and reduce the second sum to the integral representation 



^ 2/- 1 

2=0 



1 r 1 <w 







2 Jo &\y/T 




Vl + 2xd + I? 2 



2a;i? + ?9 2 

- 2) - 1 . (58) 



Evaluating the integrals in much the same way as above, 
and substituting back in (|56|l we find 



pH 



1 



1 



2V2irMl V vT 



4ttM 5 3 




2ttM| V 



(59) 



At short transverse distances x = 1 — H 2 lZ 2 /2 ignoring 
a constant and higher powers of TZ, we find the leading 
order behavior of /sd to be 



/5B(fi) 



2-kMITI 



3pH , . ft 

F3M7 



(60) 



4ttM| '"'2H- 1 '' 
The first term is precisely the 5D shock wave solution 



'5 

n is prec' 

considered in 0, [2^1 • The second term arises from the 
residual 4D graviton zero mode which persists on the nor- 
mal branch in the limit g — > 00 because the bulk volume 
remains finite. Its perturbative coupling is set by the ef- 
fective 4D Planck scale M| e ^, which can be computed 
from the normalization of the zero mode graviton wave 
function in the bulk |n 39j as 



M 4 2 e// =2M| 



dz n 3/2 = 



4Mf 
3H 



(61) 



Substituting this in the formula for the 4D shock wave 
profile at short distances (|53l) exactly reproduces the sec- 
ond term in l|60(l . Notice the relative sign difference be- 
tween the two terms on the RHS of l|60|l . This is neces- 
sary in order for both to yield an attractive force oc —V/ 
on a test particle in the shock wave background. There- 
fore the result precisely meets the expectations for 
the shock wave profile arising in a theory with both a 
5D and 4D graviton modes in flat bulk. Because the 
two gravitons are equally coupled, the theory looks effec- 
tively five-dimensional at all sub- horizon distances, where 
the inverse power of TZ in l|60|l wins over the logarithm. 
Hence indeed in this limit the effects of the 4D graviton 
are completely negligible inside the cosmological horizon, 
and the "gravity filter" has been turned off. 

This analysis clearly demonstrates how the "gravity fil- 
ter" of works even beyond perturbation theory. It is 
encoded in the behavior of the coefficients in (|29J) . which 
decrease with the index I of the Legendre polynomials. 
Because the graviton momentum on a transverse 2-sphere 
along de Sitter brane is proportional to the index of the 
polynomial, q ~ HI, the modes with q > gH ~ l/r c have 
amplitudes suppressed by q. Hence the short-distance 
divergences of the series can be no worse than a loga- 
rithm for any finite value of g. Indeed, at short distances 
TZ < t c — H 1 /^ our shock waves are exactly the same 
as in 4D GR. For very low momenta the amplitudes are 
not suppressed, and so the shock wave profile changes to- 
wards 5D behavior at large distances TZ > lf _1 /g. The 
limit g — ► 00 on the normal branch illustrates this explic- 
itly: when r c — > all modes with finite momenta remain 
unsuppressed and we recover 5D gravity. 



VI. PLANCKIAN SCATTERING 

In this section we explore the scattering of two rela- 
tivistic particles on the brane, rrik/pk *C 1, whose energy- 
momenta p\iPi reach beyond the Planck scale. In this 
regime, gravity will generically dominate over other in- 
teractions. Indeed in 4D the effective coupling of gravi- 
tons to the packets of stress-energy is a ~ Gns, where 
s = (pi — P2) 2 = AE 2 is the Mandelstam s-parameter 
and E the center-of-mass (COM) energy [U El E2 . It 
has been shown that even if it exceeds the Planck scale, 
the scattering cross section can still be reliably calcu- 
lated using the methods of quantum field theory and GR 
|40l Ell 142. I43L Classical gravitational scattering 

has been studied before by D'Eath and Pay ne . and 
the eikonal approximation calculations of [i(J, |U |43 are 
consistent with these results when the impact param- 
eter b is larger than the center-of-mass Schwarzschild 
radius ~ Gn\/~s- The problem of Planckian scattering 
and black hole formation has been revived recently in 
the context of braneworld constructions, since if the fun- 
damental scale of gravity is really low, it may be possible 
to form black holes in high energy processes in colliders 
HI El El El E3 lEa- Thus one might wonder if 
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such processes would be different in interesting ways in 
the framework of brane-induced gravity. 

However, in light of our findings above, it is clear that 
in the case of DGP with finite g, when the impact param- 
eter 6 is smaller than the scale of modification of gravity, 
6 < H^ 1 /g, Planckian gravitational scattering will be 
very well approximated by the conventional 4D GR de- 
scription. The differences are suppressed by the powers 
of m/p and H1Z, which are both very small for highly 
energetic particles at short distances. 

Let us now explicitly verify this. In the limit when 
the impact parameter is smaller than the DGP scale of 
modification of gravity, b < H^ 1 /g, the scattering is me- 
diated predominantly by the 5D modes with transverse 
momenta q > 1/6 > gH. These individual 5D modes 
couple to a brane-localized source with the effective cou- 
pling 



(62) 



as can be seen directly from <|29[1 . and is consistent with 
1]. The latter term is the "gravity filter", which cuts off 
the influence of the modes with high transverse momen- 
tum. This is in contrast to what happens in the Randall- 
Sundrum models [39| , where the modes carrying greater 
momentum compete more efficiently with the zero-mode 
graviton |35| . The couplings can be rewritten as 



(63) 



where the first factor is the usual 5D bulk coupling, with 
the factor of 2 accounting for the orbifolding, the second 
factor is the "volume dilution" ~ l/(g r c) = H, and the 
last factor is the "filter" . Hence in order for a mode me- 
diating a process on the brane with a transverse momen- 
tum q to be strongly coupled, the process should unravel 
at COM energies E where Gn e ff(q)E 2 > 1, or there- 
fore when E > ^J~^M±. This should be compared with 
the requirements for the processes in the bulk, where the 
filter does not operate, because they are mediated by de- 
localized modes, which become strongly coupled at much 
lower scales E > M 5 . 

Therefore we see that for the brane processes at en- 
ergies M^/VbH > E > M 4 , all the individual modes 
with momenta q > 1/6 are in fact weakly coupled, with 
GNeff{q) < 1- However, the resonance mimicking the 
4D graviton is strongly coupled. Indeed, this resonance 
couples to the brane-localized sources with the effective 
strength 



G 



q>l/b 



Gn eff(q) 



1 



" dq 

i/6 q 



l 



WL,7 = W HAbh 



(64) 

where the UV cutoff A is the highest transverse mo- 
mentum for which the particle scattering is meaningful, 
which we roughly estimate as A ~ M4. Thus Gn 



up to logarithmic corrections, and so at COM energies 
E > M4, indeed GnE 2 > 1, implying that the resonance- 
mediated scattering is in the strong coupling regime, 
where it should dominate over other processes much like 
the conventional 4D gravity. This description should be 
reasonable as long as the impact parameter 6 is greater 
than the Schwarzschild radius ~ GnE ~ Gn\/s ll^l, 
where black holes are expected to form [40l l42l l45l l50| , 
and where one would need to go beyond the eikonal ap- 
proximation. It is interesting however that in the regime 
of COM energies M 4 < E < M 4 /VbH, and impact pa- 
rameters 6 > GnE, which yield the COM energy inter- 
val M A < E < MiiMi/H) 1 / 3 , the strongly coupled 4D 
graviton resonance is resolved into the individual weakly 
coupled bulk modes. 

In this limit, the scattering amplitude can be extracted 
directly from the shock wave solutions 0]. For a shock 
wave profile /(f2) on the brane at z = 0, the formula 
relating the scattering amplitude in the eikonal approx- 
imation A(s,b), expressed via the eikonal ^e(s,6) — 



(65) 



m (^%^) ' an d the shock wave metric 10 is [35ll4 



P xt \ S E {s,b) 
9uu = - o(u) 

7T S 



Solving for SE(s,b) and substituting in the shock wave 
profile formula from and (|55[) . we obtain, taking into 
account our normalizations, 

5 E(s,b) = -^((l + ai H 2 6 2 + ...)ln(^) 



2H- 



1 

17? 



+ const + biHb + b 2 H 2 b z + ...J ,(66) 

where , bk are the same numerical coefficients as in eq. 
(|55|l (and 61 ~ g). Thus clearly as long as 6 <C r c = 
H^ 1 /g, the deviations away from the expression in 4D 
GR are completely negligible. Hence as we have claimed 
above, the scattering will be very well approximated by 
the analysis in the conventional 4D theory |40l l4ll l42l| . 
As a consequence one expects also that the processes of 
black hole formation in high energy collisions would be- 
have as in 4D GR |45j, |4g, |50|, controlled by the low 
energy value of the Planck scale, M4 ~ 10 19 GeV, which 
is much too high to see such processes in forthcoming ex- 
periments at the LHC. This provides further illustration 
of the efficiency of the "gravity filter" 0] . 



VII. CONCLUSIONS 

In this paper we have given a detailed derivation and 
analysis of the exact gravitational shock wave metric in 
DGP models , first presented in our earlier paper 0] . 
The shock wave solutions provide an explicit and concrete 
example of the "gravity filter" mechanism of beyond 
perturbation theory. Even on these exact solutions, the 
effects of the whole continuum of bulk graviton modes 
is dominated by a long range resonance, which mimics 
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the standard 4D graviton out to distances r c = H~ 1 /g. 
This is because the bulk modes which can transfer large 
transverse momentum on the brane q > gH ~ 1 /r c have 
amplitudes suppressed by q. Hence at short distances 
1Z < r c their influence is cut off, and the rate of diver- 
gence of the shock wave metrics can be no worse than a 
logarithm for any finite value of g, so that to the lead- 
ing order the gravitational shock waves are exactly the 
same as in 4D GR. For very low momenta the ampli- 
tudes are not suppressed, leading to IR corrections which 
change the shock wave towards the 5D form at large dis- 
tances 1Z > H~ 1 /g. We have also demonstrated that the 
leading-order Planckian scattering behaves in exactly the 
same way as in 4D GR. The differences come with the 
powers of m/p and H1Z, and are very small for highly 
energetic particles at short distances. Further, we have 
noted a spectacular new channel for energy loss into the 
bulk on the self-inflating branch, which arises when grav- 
ity is modified at exactly the brane de Sitter radius. 

We have argued earlier ^3 that our exact shock waves 
may also provide for a new arena to explore the scalar 
graviton sector in DGP. In the relativistic limit pursued 
here and in |l5j the scalar graviton has decoupled, in 



spite of the strong coupling issues that were encountered 
m naive perturbation theory 0, U El E3 

However one 

can imagine taking a very fast, but not ultra-relativistic 
observer who explores the field of a non- vanishing source 
mass. In this limit, one may treat the rest mass of the 
source as a perturbation of the shock wave geometry, 
and study how it sources the scalar graviton, and what 
happens to the graviton perturbation theory around the 
shock wave. This could reveal new aspects of modified 
gravitational physics in the IR, and help illuminate the 
status of effective 4D theory in DGP. 
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